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Abstract 

We study nucleon structure functions within the bosonized Nambu-Jona-Lasinio 
model where the nucleon emerges as a chiral soliton. We discuss the model predictions 
on the Gottfried sum rule for electron-nucleon scattering. A comparison with a low-scale 
parametrization shows that the model reproduces the gross features of the empirical struc- 
ture functions. We also compute the leading twist contributions of the polarized structure 
functions gi and g2 in this model. We compare the model predictions on these structure 
functions with data from the E143 experiment by GLAP evolving them appropriately. 

Here we discuss a possible link between two successful although seemingly unrelated pic- 
tures of baryons. On one side the quark parton model successfully describes the scaling 
behavior of the structure functions in deep inelastic scattering (DIS) processes. Scaling viola- 
tions are computed within perturbative QCD. On the other side the chiral soliton approach, 
motivated by the large Nc (number of colors) expansion of QCD, gives a fair account of many 
static baryon properties. For Nq — > oo, QCD is equivalent to an effective theory of weakly 
interacting mesons. This unknown meson theory is modeled by rebuilding the symmetry 
structure of QCD. In particular this concerns chiral symmetry and its spontaneous break- 
ing. Baryons then emerge as non-perturbative (topological) meson fields configurations, the 
so-called solitons. It is our goal to link these two pictures by computing nucleon structure 
functions within a chiral soliton model. We start from the hadronic tensor 



W, 



A e^'^<{NiP)\\rM),4H0)] \NiP)) , 



(1) 



which describes the strong interaction part of the DIS cross-section. \N(P)) refers to the 
nucleon state with momentum P and J^(0 to the hadronic current. The leading twist con- 
tribution to the structure functions is extracted from W^^(g) by assuming the Bjorken limit 



\q\ — Mjyx with |g| oo and 



-q'^/2P ■ q fixed . 



(2) 



In most soliton models the current commutator (0) is intractable. As an exception the 
Nambu and Jona-Lasinio (NJL) modelQ for the quark flavor dynamics, which can be 
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bosonized by functional integral techniques!^, contains simple current operators. Most im- 
portantly, its bosonized version contains soliton solutions]^. We confine ourselves to the key 
issues of the structure function calculation as details are given in refs.Q,^,^. 

The Nucleon from the Chiral Soliton in the NJL Model 

Here we briefly summarize the description of baryons as chiral solitons in the NJL-model, 
details are found in the review articles ||3|,|^]. We consider a chirally symmetric NJL model 
Lagrangian which contains interactions in the pseudoscalar channel. Derivatives of the quarks 
fields only appear in form of a free Dirac Lagrangian, hence the current operator is formally 
free. Upon bosonization the action may be expressed asQ 

A = Tt InA {i^ - mU^') + "^^tv (U + -2) , (3) 

where tr and Tr denote discrete flavor and functional traces, respectively. The model pa- 
rameters are the coupling constant G, the current quark mass ttiq and the UV cut-off A. 
The constituent quark mass m arises as the solution to the Schwinger-Dyson (gap) equation 
and characterizes the spontaneous breaking of chiral symmetry. A Bethe-Salpeter equation 
for the pion field {U = exp(ir • tt/Ztt)) is derived from eq. (^. Then the model parameters 
are functions of the pion mass m-j^ = 135MeV and decay constant 7,^ = 93MeV. The sole 
undetermined parameter is m. An energy functional for non-perturbative but static field 
configurations U{r) is extracted from (|3|). It is expressed as a regularized sum of all single 
quark energies e^. For the hedgehog ansatz, Uh = exp(iT • f 0(r)) the one-particle Dirac 
Hamiltonian reads 

h = a ■ p - m e-K.p{i-f5T ■ rQ{r)) , /i^f^ = e^^*^ . (4) 

The distinct level (v), bound in the background of Uh, is so-called valence quark state. 
Its explicit occupation guarantees unit baryon number. The self-consistent minimization of 
the energy functional determines the chiral angle 0(r). The so-constructed soliton does not 
carry nucleon quantum numbers. To generate them, the time dependent field configuration 
is approximated by elevating the zero modes to time dependent collective coordinates A{t) G 
SU(2) : U{r,t) = A{t)UH{'r')A^ (t). Canonical quantization of the angular velocities, Q = 
—2itr{TA^ A), introduces the spin operator J via ft = J jo? with o? being the moment of 
inertia^. The nucleon states |A^) emerge as Wigner D-functions. The action (^) is expanded 
in powers of fl corresponding to an expansion in 1/Nc- In particular the valence quark 
wave-function ^^{x) acquires a linear correction 

^,{x,t) = e--^*yl(t) I ^,{x) + J2 ^^(x)^-^j^^^ \ = e--^*^(t)Vv(^). (5) 

Here ipv{x) denotes the spatial part of the body-fixed valence quark wave-function with the 
rotational corrections included. 

^Generalizing this treatment to flavor SU(3) indeed shows that the baryons have to be quantized as half-integer 
objects. For a review on solitons in SU(3) see e.g.lql. 
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Figure. 1. The unpolarized structure functions obtained after extracting the coUective part of the 
nucleon matrix elements. Here we used m — 350MeV. 

Structure Functions in the Valence Quark Approximation 

The starting point for computing the unpolarized structure functions is the symmetric 
part of the hadronic tensor in a form suitable for localized quark fields |P, 



Snpua sgn ( 



dt / d'^xi d'^X2 e 



i{ko+qo)t 



(6) 



(2^)4 ^f^P"' 

xexp [-i{k + q) ■ {xi - X2)] {N\ [^{xi,t)tit^-f''i>{x2,0) - ^{x2, 0)t^ta"^'(a;i, t)} \N). 

Note that the quark spinors are functionals of the soliton. Here S^p^a = g^pgua + g^^gup — 
g pug pa and C = 1(2) for the structure functions associated with the vector (weak) current 
and tm is a suitable isospin matrix. The matrix element between the nucleon states (|A^)) 
is taken in the space of the collective coordinates. Eq. (|6|) is derived by assuming the free 
correlation function for the intermediate quark fields. In the limit (^) the momentum, fc, of 
the intermediate quark is highly off-shell and hence not sensitive to momenta typical for the 
soliton configuration. Thus the use of the free correlation function is a valid treatment. 

The valence quark approximation ignores the vacuum polarization in (^), e.g. the quark 
field operator ^ is substituted by the valence quark contribution (^). For m ~ 400MeV this is 
well justified since this level dominates the nucleon observables|^^. The structure function 
F2{x) is obtained from (^ by an appropriate projection]^. After computing the collective 
coordinate matrix elements all physical relevant processes are described in terms of four 
reduced structure functions /^'^ shown in figure ||. The superscript denotes the isospin of 
ti X tm while the subscript refers to forward and backward moving intermediate quarks in (|6|) . 
Although the problem is not formulated Lorentz-covariantly these structure functions are 
reasonably well localized in the interval x € [0, 1] . The contributions of the backward moving 
quarks are quite small, however, they increase with m. 

In figure we display the linear combination relevant for the Gottfried sum rule 



(Ff -Fn = -x (/|-/l)/3 



(7) 



In the Bjorken limit the Callan-Gross relation F2{x) = 2xFi{x) is satisfied. 
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Figure. 3. The polarized structure functions gi and 52 after projection and QCD evolution. Left 
panel: The dashed (dotted) line denotes the (projected) low scale model prediction. 



and compare it to the low-scale parametrization of the empirical data|lC]. This is obtained 
from a next -to-leading order QCD evolution of the experimental to a low-energy regime typ- 
ical for soliton models. The model reproduces the gross features of the low-scale parametriza- 
tion. Moreover the integral of the Gottfried sum rule 



sg = r^{F^^-Fn = \ 

Jo X { 



0.29 , m = 400MeV 
0.27 , m = 450MeV 



(8) 



reasonably accounts for the empirical value Sq = 0.235±0.026[11]. In particular the deviation 
from the nai've value (1/3)[12| is in the direction demanded by experiment. In figure ^ we 
also compare the model prediction for the polarized structure function gi{x) to its low-scale 
parametrization pO[] . Apparently they agree the better the smaller the constituent quark mass. 
This behavior is opposite to the unpolarized case. No low-scale parametrization is available 
for the polarized structure function g2{x). Therefore we have projected the corresponding 
prediction onto the interval x € [0, 1] ||l3|,|l^ and subsequently performed a leading order 
QCD evolution to the scale of the experiment B]. Figure shows that the resulting polarized 



structure functions reproduce the empirical data quite well, although the latter have sizable 
errors. 

Note that for consistency with the Adler sum rule also the moment of inertia must be 
restricted to the valence quark contribution For m = 350MeV this fortunately is almost 
90%. Analogously the Bjorken, Burkhardt-Cottingham as well as the axial singlet charge 
sum rules are confirmed within this model treatment |^] . 

Conclusions 

We have presented a calculation of nucleon structure functions within a chiral soliton 
model. We have argued that the soliton approach to the bosonized version of the NJL- 
model is most suitable since (formally) the current operator is identical to the one in a free 
Dirac theory. Hence there is no need to approximate the current operator by e.g. performing a 
gradient expansion. Although the calculation contains a few (well-motivated) approximations 
it reproduces the gross features of the empirical structure functions at low energy scales. This 
is true for the polarized as well as the unpolarized structure functions. Future projects will 
include the extension of the valence quark approximation, improvements on the projection 
issue and the extension to flavor SU(3). 
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